Quantum trajectory analysis for electrical detection of single-electron spin resonance 
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A Monte Carlo simulation on the basis of quantum trajectory approach is carried out for 
the measurement dynamics of a single electron spin resonance. The measured electron, which 
is confined in either a quantum dot or a defect trap, is tunnel-coupled to a side reservoir and 
continuously monitored by a mesoscopic detector. The simulation not only recovers the observed 
telegraphic signal of detector current, but also predicts unique features in the output power 
spectrum which are associated with electron dynamics in different regimes. 
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PACS numbers: 03.67.Lx, 85.35.-p, 73.63.-b, 76.30, 



In recent years, largely being stimulated by the in- 
terest of solid-state quantum computation, measurement 
of single electron spin has become an intensive research 
subject. The main difficulty of spin measurement lies 
in the inherent weakness of magnetic interaction, mak- 
ing its direct detection quite challenging. The current 
state-of-the-art techniques include, for instances, scan- 
ning tunneling microscopy 0, magnetic resonance force 
microscopy (MRFM) , sub-micron Hall magnetometer 
Q, and spin-to-charge conversion Q, etc. Experimen- 
tally, detection of single electron spin has been illustrated 
by optical means ElE!' MRFM technique 0, and elec- 
trical methods Ulij. In particular, the idea of electri- 
cal approach is based on the spin-to-charge conversion, 
which maps the spin-states onto detectable charge states. 
In Refs. @ and0 single shot read-out of the charge state 
is demonstrated in terms of single realization of contin- 
uous measurement. In this work, we present a Monte 
Carlo simulation for the dynamics of continuous electri- 
cal measurement, and calculate the output noise spec- 
trum. This study provides a fundamental description 
for the stochastic nature of state evolution, which stems 
from the interplay of the manipulation of coherent field 
and the (simultaneous) influence of incoherent reservoir. 
Also, the predicted unique features in the noise spectrum 
manifest distinct electron dynamics in different regimes. 

In close connection with the experiments [3,13 > we con ~ 
sider a model system of quantum dot in the Coulomb 
blockage regime, which is tunnel-coupled to a side elec- 
tronic reservoir and continuously monitored by an elec- 
trical detector, say, the quantum point contact (QPC), as 
schematically shown in Fig. 1. The single-electron level 
in the dot is split by an external magnetic field Bq, i.e., 
with the Zeeman energy of ej_ — ej = g/is-Bo = A, where g 
is the Lande-g factor of electron spin and is the Bohr 
magneton. The spin-up and spin-down states are coupled 
by an oscillating magnetic field -Besr — -Bi cos(wt), ap- 
plied perpendicularly to the field Bo, which leads to the 
well known electron spin resonance (ESR) with Rabi fre- 
quency wr = gfisBi. In the Coulomb blockade regime, 




FIG. 1: Schematic diagram of a quantum dot being tunnel- 
coupled to a side reservoir and measured by a quantum point 
contact (QPC). The single-electron level in the dot is split 
by an external magnetic field and the two spin states are 
coupled by an electron spin resonance (ESR) magnetic field. 
The tunneling of single electron into and out of the dot would 
alter the transport current through the QPC detector. 



the dot is occupied by at most one electron. Thus, the dot 
electron Hamiltonian in the presence of magnetic fields 
reads i/dot = f (d\d i ~d\d ] ) + ^f(d\d ] e- wt + d\d ie lult ), 
in which the first part describes the Zeeman splitting, and 
the second one represents the electron spin resonance. 
Following the QPC measurement model proposed by 
Gurvitz |lf)l | , the entire measurement setup shown in Fig. 
1 is described by the Hamiltonian H = Ho + Ht + -Hint , 
with 



Ho — Hdot + H ICS + Hpc 
H T = J2[t q cl a d a + H.c] 



qa 



H n t = ^(Tqk + Xqkno)a Lk a Rq + H, 

k.q 



(la) 
(lb) 

(1c) 



Here the Hamiltonian H ICS = t q c^ qcy c qa stands for 
the side electronic reservoir, and Ht describes its tunnel- 
coupling to the quantum dot. a =\,\ denote the two 
possible orientations of electron spin. For the QPC de- 
tector, H PC = J2k e k a lk a Lk + J2q e q a R q a R<] stands for 
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the two (left and right) reservoirs, and H lnt describes tun- 
neling between them. Notice that the tunneling ampli- 
tude explicitly depends on the electron occupation of the 
quantum dot, which is modelled by T q k + Xqk^o, where 
no = 1— d\d a is the vacancy operator of the quantum 
dot. 

Reduced description for the measurement process. - 
A significant aspect of quantum measurement is to in- 
vestigate the back-action of measurement device on the 
measured system, which results in the dephasing and 
relaxation. Technically, this can be realized by tracing 
out the microscopic degrees of freedom of the measure- 
ment device. Physically, this reduced description cor- 
responds to ensemble results, i.e., all the measurement 
records are averaged. For present setup, the side reser- 
voir and the detector can be treated as two independent 
environments acting on the (measured) quantum dot. 
As a model description, we assume the tunneling am- 
plitude Tfe g to be real and independent of the reservoir- 
state "/eg", i.e., H[ nt = QF, where Q = T + x n o an d 
F = J2kq( a \,k a Rq + H.c). For the coupling between 
the dot and the side reservoir, we re-express the cou- 
pling Hamiltonian as Ht = ^2 a {d a fl + H.c), where 
f\ = J2 q tq c \a- By treating Ht and Hi nt as perturbation 
and on the basis of the second-order cummulant expan- 
sion, an equation for the reduced density matrix can be 
derived as [UGl 



i£op- -AQ,Qp- pQ + ] 

^{[dlD^p-pD^+H.c] 



(2) 



Here, the first term C$p = [H^ t,p] describes the coher- 
ent manipulation of the dot electron (by the ESR mag- 
netic field), while the second and third terms stem from 
the back-action of the detector and the random tunnel- 
ing caused by the side reservoir. In the back-action term, 
the operator Q = Q(+) + Q ( ~\ and Q(±) = C {± \C )Q, 
where the superoperators (Co) are the Fourier trans- 
form of the correlation functions C?W(t) = (F(t)F(0)) 
andC ( -)(i) = (F(Q)F(t)). Under the wide-band approx- 
imation, C^(Co) can be explicitly carried out as [Tl| : 

c (±) (£o) = v 



. Here, -q = 2itgig ri 

=-(£o±eV) 

with gi and g r the energy-independent density of states 
for the two reservoirs of the QPC detector. f3 = l/(ksT) 
is the inverse temperature, and eV = pl — P>r i s the ap- 
plied voltage across the detector. Similarly, for the side 
reservoir, we obtain \V2}. D { a ±} = T a ni ±] ' (e a )d a . Here 

ni + \ecr) = n a (e a ) and n£ '(e CT ) = 1 - n a (e a ), with n a 
the Fermi distribution function, = 2irg c \t q \ 2 , in which 
g c is the density of states of the side reservoir electrons 
at the energy e a . 

At the limit of large bias voltage , Eq. (J2J can be 

further simplified to the Lindblad-type master equation 



Here we have also restricted our study to zero tempera- 
ture. The superoperator T> is defined as D[r]p = J[r]p — 
A[r\p, where J\r\p = rpr\ and *4[r]p = ^(r^rp + prh*). 
The jump operators Si = |0)(| | and S 2 = \ T)(0| de- 
scribe the tunneling between the dot and the side reser- 
voir. The homodyne-type jump operator D[T + ^no], 
which describes the back-action of the detector, is asso- 
ciated with the tunneling amplitudes through the QPC 
with and without an electron in the quantum dot, by the 
corresponding tunneling rates |T| 2 = 7]\T\ 2 eV = D, and 
\T + X\ 2 = rj\T + x?eV = D'. Note that, differing from 
the setup of a charge qubit measured by a QPC [l(|, here 
the QPC's back-action has minor effect on the measured- 
electron dynamics. The reason is that, in the occupation 
state basis of the quantum dot, D[T + Xno]p only con- 
tributes to the off-diagonal elements, which are decoupled 
from the diagonal ones. That is, the back-action only 
causes dephasing between the "occupied" and "empty" 
dot states. But no phase coherence is introduced initially, 
and neither will it be generated in the later on evolution. 
Thus the QPC has no influence on the electron dynamics 
under study [Tij . 

Straightforwardly, the (ensemble) measurement cur- 
rent through the detector can be calculated via 



I = IoPqo + hpn, 



(4) 



P = - ^A ) p + ^ i 2?[s 1 ] i0 + ^ T p[ 1 s 2 ]p 

+ V[T + Xn }p = £p(t). 



(3) 



where Iq = eD' and I\ — eD are, respectively, the cur- 
rents through the QPC in the absence and presence of 
an electron in the quantum dot. Note that Eq. is pre- 
cisely the result presented in Ref. ^[ In the following 
we will see that it is indeed an ensemble average of the 
continuously measured currents from a large number of 
individual realizations. 

Quantum trajectory description. — The above reduced 
density matrix description corresponds to average over all 
the measurement records of many individual realizations. 
However, in most cases the original readout data of con- 
tinuous measurement of a single quantum system is the 
record from a single measurement realization. For in- 
stance, the experimental result of Ref. reveals random 
telegraph signal (RTS), with detector current jumping 
between two discrete values stochastically. 

Proper description for this single measurement realiza- 
tion is the quan tum trajectory theory developed in quan- 
tum optics [Tfl [13] , which describes the system evolution 
conditioned on the information continuously acquired by 
the detector. The central ingredient of this formalism 
is the conditional master equation (CME), which is an 
unravelling of the unconditional master equation Eq. 

For the tunneling events between the quantum dot and 
the side reservoir, two stochastic point variables dN\{t) 
and dN2(t) (with values either or 1) are employed to 
denote, respectively, the electron numbers that tunnel 
from the dot to the side reservoir and vice versa, during 
the infinitesimal time interval dt. While for the QPC, 
at the limit of zero temperature and large bias voltage, 
the electron only tunnels from the source (left reservoir) 
to the drain (right reservoir), which is denoted by an 
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alternative stochastic variable dN c (t). Following Ref. IT3I 
we obtain 



dp c {t) 



dt{-i£ - A[T + Xn ] 
-T^S,] - T^A[S 2 ] 
+V c {t) +Vx(t) + V 2 {t)} p c (t) 



+dN 1 (t) 
+dN 2 {t) 
+dN c (t) 



1 



T^J[S 2 ] 

V 2 {t) 
J[T+Xn ] 
■Pc(t) 



Pc(t) 
Pc(t) 



Po(t), (5) 



in which the stochastic variables dNi(t) (i = c, 1,2) 
have the ensemble property E[dNi(t)] = Vi(t)dt = 
Ti{J[ri]p c (t)}dt, with r c = T + Xuq, r\ = w / TjS'i and 



r 2 = VM^2, respectively. Here £[-••] denotes an en- 
semble average of a large number of stochastic processes. 

Within the quantum trajectory approach, the current 
through the QPC is given by 



I(t) =edN c (t)/dt. 



(6) 



Here we show that its ensemble average recovers the cur- 
rent given by Eq. Q. With the help of the CME Eq. JSJ), 
the stationary-state current can be evaluated as 



Ioc = eE[dN c (t)/dt] t ^oo 

= eTi {J[T + Xn ]p(oo)} 

D + (D' - D) 



r 2 



3u% 



AS 2 



(7) 



where p(oo) is the ensemble stationary state calculated 
by Eq. 0. Here, we have assumed — T[ = T. We 
notice that Eq. (JZJ) is precisely the same as derived in 
Ref. [3 based on Eq. (gj. 

For the typical mesoscopic detector QPC, the condi- 
tion D ^ \D' — D\ implies a large background cur- 
rent. Accordingly, it is desirable to convert the above 
point process description into a diffusive one fl3| . First, 
the terms in Eq. ©, {~A[T + Xn ] + T c (t)}p c (t)dt + 

dN c (t) \ J %+u? o] ~ lj Pc(t), are replaced by V[T + 
Xn Q } Pc (t)dt + X[n oPc (t) + p c (t)n Q - 2(n )p c (t)]£(t)dt, 
where t;(t) is a Gaussian white noise characterized by 
E[£(t)] = and E[£(t')€(t)] = S(t - t'). Then, the detec- 
tor current I(t) of Eq. 10 in the diffusive regime becomes 



I(t) - I = AI[p c 00 (t) - Pii(t)}/2 + y/Sofim, (8) 

where I = (Iq + Ii)/2 , AI = Iq — Ii, and the background 
shot noise So = 2e/oo. Here we introduced the notation 
Pij(t) = (i\p c \j), with u (i,j) = (0,1)" characterizing the 
number of electron occupied in the dot. 

From Eq. (JSJ we find that the detector current I(t) 
basically follows the occupation probabilities and 
contains a noise term £ (t) that purely arises from the in- 
trinsic fluctuations in the detector, being independent of 
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FIG. 2: Detector current from a single particular realization 
of quantum measurement, where the feature of random tele- 
graph signal (RTS) which jumps between two discrete current 
values is evident. In the simulation, the resonance condition 
(i.e. 8 — A — uo — 0) and weak coupling between the dot and 
the side reservoir (by setting V — Q.Iujr) are assumed. The 
simulated random current I(t) in (a) is filtered by using a rect- 
angular window with width r = 0.9/ujr, and in (b) is further 
averaged over the detector's background current fluctuations. 



the measured system. Particularly, in the Markov limit, 
the detector noise subjects to Gaussian distribution as 
characterized above. For noise strength weaker than the 
signal, single particular realization of measurement will 
be evident. Since the fluctuation grows when I(t) is ex- 
amined at smaller time scales, proper averaging over time 
("low-pass filtering") is necessary [l8|. Based on Eq. JSJl, 
we thus average I(t) by using a rectangular filter window 
with time width r, I(t) = r^ 1 f*_ T I(t')dt' . The width 
of the filtering window should be larger than the average 
time interval of two successive electron tunneling events 
through the QPC junction, and smaller than the average 
occupation/empty time of the quantum dot. This makes 
the output current not too noisy in order to extract the 
information of the measured system. The numerical re- 
sult is shown in Fig. 2 (a), where the width of the filtering 
window is chosen as r = 0.9/Vr. As a manifestation of 
the telegraph signal, the current stochastically switches 
between two discrete values. The low and high current 
plateaus correspond to the quantum dot with and with- 
out an extra electron in it. The average widths of the 
plateaus are determined by the average occupation prob- 
abilities Pii(t). Also, the detector shot noise would cause 
current fluctuations around the plateaus. If the detec- 
tor's background current fluctuations are further aver- 
aged, the telegraph signal will be more evident as shown 
in Fig. 2 (b), which agrees well with the experimental 
observation reported in Ref. 0, in the regime T <C u>r. 

As we have seen in Fig. 2, the single realization of 
the output current largely reflects the instantaneous oc- 
cupation of the quantum dot, which is determined by 
the interplay of the coherent ESR driving and the in- 
coherent tunneling. Specifically, we can easily imagine: 
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FIG. 3: Stationary noise power spectrum of the detector cur- 
rent. (a)-(d) correspond to different coupling strengths F be- 
tween the quantum dot and the side reservoir. 

(i) in the weak tunneling regime, the dot is either occu- 
pied or unoccupied for a long time on average; (ii) in the 
intermediate tunneling regime, the dot is occupied for 
longer time than unoccupied; and (iii) in the strong tun- 
neling regime, the dot is largely occupied. Accordingly, 
the output current would follow these occupations and 
exhibit the corresponding behaviors. Rather than show- 
ing the entire landscape of the dot-electron dynamics by 
the output current, in the following we show that it can 
be elegantly manifested in the output power spectrum. 

Noise spectrum. — Following Refs.ll9l and 2& as deriv- 
ing the ensemble-averaged current Eq. J7J, starting with 
the current formula Eq. JfJJ which is defined in the "point 
process" regime can also lead to the current correlation 
function in the stationary state as 

G(r) = {E[J(ttr)I(i)]-B[J(tt«i)]U 
= e/oo^r) +e 2 (I/ - D) 2 

x {Tr [n e CT J[n Q ] Poc ] - Tr[n Poo] 2 } . (9) 

It should be noted that although this expression was de- 
rived in the "point process" regime, it is still applicable 
in the diffusive regime, since it is an ensemble averaged 
result. By noting G(— t) = G(t), the stationary noise 
power spectrum can be computed by 

/oo 
dTG(r)e- lWT . (10) 
-oo 

Numerical results of the noise spectrum are shown in Fig. 
3, where the dot-electron dynamics under the interplay 
of the coherent ESR driving and the incoherent tunnel- 
ing is clearly manifested. In the weak tunneling regime 
(i.e. r <C u>r), compared to the electron spin-up and 
spin-down flipping time in the dot, the electron tunnels 
into and out of the dot slowly. As a consequence, the dot 
occupation only changes occasionally with long-time sep- 
aration, which results in the dominant sharp Lorcntzian 
peak around the zero frequency, see Fig. 3(a). This fea- 
ture indicates nothing but the telegraph noise. With the 
increase of T, the in-dot coherent flipping will eventually 



FIG. 4: ESR-induced coherent oscillation of the electron oc- 
cupation probability in the quantum dot (solid curve). The 
non-diagonal element pn describes the phase coherence of the 
ESR (dashed curve), where the decay is caused by the tunnel- 
coupling with the side reservoir (but not the back-action of 
the QPC detector). The adopted parameters are the same as 
for Fig. 3(b). 



play important role on the electron's dynamics, which 
leads to the gradual suppression of the central sharp peak 
(i.e. the telegraph noise spectrum), and makes the co- 
herent peak around the Rabi frequency more evident, as 
shown by Fig. 3(b) and (c). Finally, as T ^ gjr, the inco- 
herent tunneling between the dot and the side reservoir 
takes place so fast that the dot is almost occupied for all 
time, which greatly reduces the excess noise as shown in 
Fig. 3(d). The considerably flattened peak shown there 
simply reflects the sudden change of the dot occupation, 
which resembles the white noise process. 

The above behaviors of the noise spectrum are closely 
related to the output current characteristics in different 
regimes, as discussed previously in the context of Fig. 2. 
Below we further discuss the coherent peaks in Fig. 3: (i) 
We would like to point out that the parameters D and 
D' are related to the background noise Sq = 2e/oo (the 
first term of Eq. (9)), meanwhile D — D' would scale 
the magnitude of the excess noise which is caused by the 
electron dynamics in the quantum dot (see the second 
term of Eq. (9)). However, D and D' have no effect on 
the intrinsic structure of the noise spectrum, (ii) The 
height of the coherent peak in the noise spectrum is not 
very sensitive to the tunneling rate T, but its evolution 
from "birth" to "death" by increasing T is evident. In- 
stead of its absolute value, the significant feature is the 
relative magnitude of the central (zero frequency) peak 
and the side (coherent) peaks, which essentially reflects 
the underlying dynamics, say, being incoherent or par- 
tially coherent, (iii) The position of the coherent peak is 
around the Rabi frequency ujr. This is understandable, 
because it is the Rabi oscillation which causes the "coher- 
ent" oscillating behavior of the occupation probability. 
As we have pointed out, no quantum coherence is intro- 
duced and caused between the occupied and unoccupied 
dot-states. However, the coherent flipping of the spin-up 
and spin-down states of the electron inside the quantum 
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dot would induce the coherence- like oscillation of the oc- 
cupation probability, with the same Rabi frequency, as 
shown in Fig. 4. The QPC can probe this short-time os- 
cillating behavior and results in the "coherent" peaks in 
the output noise spectrum. Precisely speaking, the peak 
position does not locate exactly at the Rabi frequency. 
This can also be understood from the analytical result 
of the noise spectrum in the measurement of solid-state 
charge qubit by the QPC |21|. However, the slight shift 
of the peak position to higher frequency with T does not 
provide any interesting information. 

In summary, we have presented a Monte Carlo simu- 
lation for the single realization of continuous detection 
of single electron spin. The measurement scheme, say, 
based on the spin-to-charge conversion and performing 
the electrical detection by a mesoscopic detector, is of 



high interest for quantum computing. The present study 
provides a theoretical description for the measured RTS 
result 0, E| , which may represent an experimental break- 
through of single electron spin detection in solid-states. 
The predicted feature in the output noise spectrum re- 
veals an interesting interplay of coherent oscillation and 
incoherent tunneling. Its possible observation in experi- 
ment has the value of directly illustrating the single elec- 
tron quantum dynamics in different regimes, i.e., from 
(partial) quantum to classical. 
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